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In an effort to investigate the link between failure mechanisms and the geometry of fractures of 
compacted grains materials, a detailed statistical analysis of the surfaces of fractured Fontainebleau 
sandstones has been achieved. The roughness of samples of different widths W is shown to be self 
affine with an exponent C, = 0.46 ± 0.05 over a range of length scales ranging from the grain size d 
up to an upper cut-off length £ ~ 0.15V4 7 . This low ( value is in agreement with measurements on 
other sandstones and on sintered materials. The probability distributions nsz(Sh) of the variations 
of height over different distances Sz > d can be collapsed onto a single Gaussian distribution with 
a suitable normalisation and do not display multifractal features. The roughness amplitude, as 
characterized by the height— height correlation over fixed distances Sz, does not depend on the 
sample width, implying that no anomalous scaling of the type reported for other materials is present. 
It is suggested, in agreement with recent theoretical work, to explain these results by the occurence 
of brittle fracture (instead of damage failure in materials displaying a higher value of ( ~ 0.8). 

PACS numbers: 62.20.Mk,46.50.+a,68.35.Ct 



I. INTRODUCTION 

The characterization of the geometry of fracture sur- 
faces is of considerable interest for a broad range of prac- 
tical applications ranging from the hydrology of fractured 
reservoirs [3, Q to the derivation of friction laws 0] . Their 
statistical analysis is also relevant to the understanding 
of the physics of fracture: The roughness may indeed be 
expected to reveal failure mechanisms occurring at the 
microstructure scale Q. It is now well established that 
fracture surfaces satisfy a scaling invariance known as 
self-affinity 0, H, 0, 0, 0, E3 ■ A self-affine profile is one 
that is statistically invariant under the scaling relation: 
z — > \z and h(z) — > X^h(z) where the z axis is in the 
mean fracture plane and normal to the direction of crack 
propagation, h(z) is the surface height profile and £ is the 
self affine exponent. The scaling invariance implies that 
the probability ns z (5h) to observe a variation of height 
5h over an horizontal distance 8z follow the scaling law: 

X C Tr XSz (X C Sh) = irsz(6h) (1) 

The prefactor A 1 " results from the normalization of 
irs z (Sh). As a direct consequence, the height— height cor- 
relation function Ah(8z) satisfies: 

Ah{Sz) = ((h(z + Sz) - hiz)) 2 ) 1 / 2 oc 6z c (2) 
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In many materials including rocks like granite and 
basalt 0, H, Q , the roughness exponent £ was found to 
be close to 0.8, suggesting at first that it was a universal 
value 5, 6]. However, different exponents were reported 
for Berea sandstone £ ~ 0.4 — 0.5 and then on syn- 
thetic sandstone made of sintered glass beads [l2T| , In 
this latter case, the exponent £ is equal to 0.40 ± 0.04 
independent of the porosity in the range investigated 
(3 % < (p < 26 %). This difference raises important 
fundamental questions regarding the physics of fracture, 
namely whether £ is determined by the structure of the 
material or by the failure mechanism. 

In support of the second possibility, a recent work [l3[ 
suggests that the value of £ depends on the existence of 
a process zone where damage processes take place. On 
the one hand, at length scales smaller than the size of 
the process zone, fracture surfaces would develop from 
the growth and coalescence of microdefaults (cracks or 
voids): This damage failure (also called quasi-brittle fail- 
ure) would lead to an exponent ( ~ 0.8. On the other 
hand, at length scales larger than the process zone, the 
fracture surfaces would result from the continuous prop- 
agation of a crack. For this brittle failure, the exponent 
would be ( ~ 0.4, reflecting, at the continuum scale, the 
toughness distribution of the microstructure scale. Two- 
dimensional numerical simulations realized under these 
two types of assumptions confirm these results by pre- 
dicting lower £ values for brittle fractures [l4[ than when 
damage is introduced fTsT ]. 

While most studies dedicated to the roughness of 
cracks dealt with surfaces with £ ~ 0.8 [1, |g, 0, H, H, 
Hoi . [H, [I?], [H, EH, the present work investigates these 
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problems through an extensive analysis of the statistical 
properties of the roughness of fractured Fontainebleau 
sandstone for which £ ~ 0.4 — 0.5. 

In addition to the self-affine exponent £> another im- 
portant statistical characteristic is the n distribution 
of the surface height variations: On granite surfaces 
(C = 0.8), it was found to be Gaussian [19| and a single 
self-affine exponent was needed to describe the correla- 
tion of the surface heights along the crack front direction. 

Another important point is the dependence of the 
roughness statistics on the size of the samples and on 
the distance to the initiation zone. In the transient 
propagation region of fractures initiated from a straight 
notch, the amplitude of the roughness is observed to de- 
pend at all length scales on the distance to the initia- 
tion [lj| [I?], ■ This variation may be described by a so 
called anomalous scaling law involving a new exponent 
Cgiobai > C referred as the global roughness exponent. A 
consequence of anomalous scaling is expected in regions 
far from the initiation and where the roughness is fully 
developed [13, The roughness amplitude in a win- 
dow of size Sz normal to the crack propagation would 
then scale with the width W of the fracture surface as: 

Ah(Sz) oc 6z <: W c <' l ° bal ~ t (3) 

In the present study, a statistical analysis of the same 
type is performed on fractured Fontainebleau sandstone 
samples of various widths W (perpendicular to the direc- 
tion of crack propagation). For each surface, the rough- 
ness exponent £ and the range of length scales over which 
a self-affine description is valid are inferred from the vari- 
ations of the correlation function Ah with Sz. Their de- 
pendence on the width W of the samples and on the lo- 
cation on the fracture surface is investigated in order to 
look for the anomalous scaling features observed on other 
materials. The probability distribution of the variations 
of height over given distances Sz are also investigated, in 
order to evaluate the magnitude of multiscaling effects 
and their dependence on Sz. 

The differences between these results and the charac- 
teristics reported previously for other materials will be 
discussed and we shall seek to relate them to the failure 
mechanisms taking place at scales larger than the grain 
size in Fontainebleau sandstone. 



II. SAMPLE CHARACTERIZATION AND 
EXPERIMENTAL SETUP 

The X-Ray analysis of the Fontainebleau sandstone 
samples used in the present work indicated that they 
contain 99 % quartz. The size of the grains range from 
100 fim to 500 fixsx with a mean value d ~ 270 fim. 16 
SEM micrographs covering a total area of 20 mm x 20 mm 
were used to measure the porosity. The pictures are 
thresholded at a gray level adjusted so as to discrimi- 
nate between grains and voids. The porosity, defined as 
the ratio between the void and total areas, is found to be 
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FIG. 1: Schematic view of the toroidal sample (outside radius 
60 mm, inside radius 25 mm) and of forces applied during the 
modified Brazilian test. 



10 ± 1 %. 

Slices of widths ranging from W — 7.8 to 51.5 mm were 
sawed off from Fontainebleau paving stones. Toroidal 
samples were obtained by using hole saws to core out the 
slices (see Fig. [1]). This was achieved by first piercing 
a circular hole with a 25 mm diameter saw and then a 
second hole of diameter 60 mm using the same spinning 
axis. Four slices of widths W = 7.8, 14.1, 26 and 51.5 mm 
were machined and placed in a hydraulic press to per- 
form modified Brazilian fracture tests [2l|. The latter 
is achieved by applying a uniaxial stress at two opposite 
points on the outside of the torus. The load is increased 
until two symmetrical cracks propagate from the central 
hole toward the outside where the compressive forces are 
applied. 

After the failure, one of the fractured blocks is secured 
to a computer controlled horizontal displacement table. 
The samples are moved stepwise using steps of lengths 
Az (See Tab. U for a detailed list of the actual values.) 
At each new position, a sensing needle is displaced ver- 
tically by a computer controlled microstep motor with 
a resolution of ~ 1 /im. The contact is detected by a 
vertical deviation of the needle. The repeatability of the 
measurement, verified by scanning several times the same 
profile, is 1 /im. After the detection of the contact, the 
sensor is raised by 200 /im before the sample is moved lat- 
erally again. The sequence is repeated n z times. Then, 
the sample is displaced in the perpendicular direction 
along the a;-axis by a step of length Ax and the scanning 
is repeated again. Finally, one obtains surface maps in- 
cluding n z x n x points. Figure [5] shows one of these maps 
that displays troughs and bumps with, for the largest, a 
typical amplitude of the order of one millimeter. In the 
following, we discuss a detailed statistical analysis of the 
spatial correlations of the roughness. 
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FIG. 2: Topographic image (248 x 252 pixels) of a fractured 
sandstone surface with a crack propagation along the x-axis. 
Sample width: W — 26 mm (along the z-axis). 
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TABLE I: Characteristic parameters of fracture surface geom- 
etry. W (mm): specimen width. Az, Ax (fim): measurement 
steps along the z and x. n z and n x : number of recorded points 
along the z and x directions. £ and £ft: self-affine exponents 
of profiles oriented in the z direction measured respectively 
using the correlation function and the Fourier power spec- 
trum. 



profiles are self-affine and characterized here by a rough- 
ness exponent £ = 0.43. 

The robustness of this self-affine description was verified 
by using other statistical methods such as the Fourier- 
analysis of the profiles [22|. The exponents obtained in 
this way are given in Tab. U and are globally indepen- 
dent of the method: Regarding the dependence on W, a 
global increasing trend is observed but the total variation 
is very weak, less than 0.1, which is of the same order as 
the interval of confidence on the value of one exponent 
estimated for a given width. The self-affine exponent ob- 
tained after averaging over all data is £ = 0.46 ± 0.05. 
Applying the same analysis to profiles parallel to the x 
direction, i.e. the direction of crack propagation, leads 
to a self-affine exponent equal to 0.49 ± 0.05. These ob- 
servations are consistent with measurements performed 
on other natural and artificial sandstones [ll|, [12, 13 • 

For Sz larger than £, the variations of Aft level off in- 
dicating that the surface is flat at large scales. In order 
to test the robustness of the value of £ with respect to 
possible bias introduced by the statistical analysis, the 
computation was repeated on data sets including only 
one out of two data points. The value of the cut-off 
length remained the same. Finally, the same procedure 
was applied in different regions of the surface. Each of 
the surfaces was divided into three regions of equal width 
W/3. Two of these subsurfaces are thus on the sides 
of the sample while the third one covers the central re- 
gion. In Fig. [31 correlation functions Ah computed on 
the central (x) and on one of the lateral (■) subregions 
are overlaid. All curves fall on top of each other, showing 
that the crossover length £ is independent of the region 



III. EXPERIMENTAL RESULTS 
A. Characterization of the surface roughness 

We first analyze profiles parallel to the z-axis, i.e. nor- 
mal to the direction of crack propagation, and located far 
enough from the initiation so that the roughness proper- 
ties are statistically stationary. The scaling properties of 
these profiles are characterized by their ID height— height 
correlation function Ah(Sz) defined by Eq. (J2]). 

Figure [3] represents in a log-log scale the variations of 
Ah as a function of the distance Sz for one of the sam- 
ples. Two distinct behaviors are observed: For Sz < £, 
Aft, varies as a power law of Sz (straight line in a log-log 
representation) while, for Sz > £, Aft remains roughly 
constant ("plateau domain"). The crossover length scale 
£ is, here, defined as the abscissa of the intersection be- 
tween the power law fit (solid line in Fig. [3]) and the 
plateau variation (dotted horizontal line) . The local slope 
of the correlation function (see inset of Fig. [3J is constant 
except for Sz values smaller than ~ 100 /im: This corre- 
sponds roughly to the grain size d which represents there- 
fore a lower boundary of the self affine domain. These 
results indicate that, at intermediate length scales, the 




Sz (/vm) 

FIG. 3: Log-log representation of the height— height correla- 
tion function Aft as a function of Sz computed along the crack 
front direction (z-axis) for a sample of width W = 7.8 mm. 
Straight line: Linear regression of the data for Sz < £ with 
slope £ = 0.43. Inset: Local slope of the correlation function. 
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of the fracture surface. Therefore, £ is a relevant length 
scale of the problem which reflects neither an influence 
of the sides of the sample nor of the sampling period of 
the surface maps. 

Before studying the scaling behavior of £ with the sam- 
ple width, we analyze now the distribution tt$ z of the 
variations of height. 



B. Statistical distribution of the fluctuations of 
height 

For each sample, a scan including at least 30, 000 mea- 
surement points has been performed. Such a large data 
set allows one to compute with a good accuracy the prob- 
ability distributions -k$ z (Sh) for values of 5z ranging from 
the typical grain size up to several times the grain diam- 
eter. Fig. 2] displays such probability distributions cor- 
responding to several Sz values ranging from 300 /jm to 
1.3 mm, i.e. larger than the mean grain diameter d. 

For self-affme profiles, one expects n$ z to satisfy 
Eq. ([1]), valid for any scaling parameter A. Using the 
particular value A = 1/Sz, this equation becomes : 



0.5 



TT6z{()h) = ( — | 7T 



1 



C 



Sh 



(4) 



Therefore, using the normalized variable u = 6h/Sz^, 
all the distributions Sz^t^sz should collapse onto a same 
master curve 7r(u), at least for Sz values pertaining to 
the self-affine regime. Such a collapse is indeed observed 
in Figure |U The same data as in the insert have been 
plotted in these normalized coordinates using the rough- 
ness exponent £ = 0.43 measured from the scaling of the 
correlation function (see Fig. [3}. 

One of the key consequence of this precise collapse is 
that all moments ((h(z + Sz) — h(z)) p ) 1 J p of order p > 1 
of the irsz distributions must scale as Sz^ . This is a clear 
evidence of the lack of multiscaling of the surface inves- 
tigated. In other words, the single exponent £ is enough 
to describe the full statistics of the variations of height 
on the fracture surface along the z-axis, at least in the 
range of length scales d < Sz < £. 

Let us now focus on the master curve n(u) which has 
the characteristic bell-like shape of Gaussian distribu- 
tions. This is more precisely demonstrated in the semi- 
logarithmic plot of Fig. O in which a Gaussian curve (con- 
tinous line) with a second moment equal to one is ob- 
served to coincide with the experimental data sets cor- 
responding to Sz > 300/im, i.e. Sz > d. The statistical 
properties of the surface are therefore fully described by 
this Gaussian distribution and by the single roughness 
exponent £ = 0.43. 

We have also represented the ^-distributions corre- 
sponding to length scales smaller than the mean grain 
diameter d — 270 /im, namely Sz = 75 and 150 /im. 
They clearly do not collapse with the other curves al- 
though their second moment was shown to scale roughly 
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FIG. 4: Probability distribution of height variations over a 
distance Sz normalized using Eq.Q with £ = 0.43 as ob- 
tained by the analysis of the second moment Ah of the height 
distribution (See Sec. MI Af) . 5h and Sz have been normal- 
ized by a same length scale length scale £ = 35 /im chosen 
so that the second moment of the normalized distributions is 
equal to 1. Inset : Plot in linear coordinates of ns z (5h) for 
300 /im < Sz < 1.3 mm. 



as Sz'' in this range of scales (see Fig. [3]). It means that, 
in the domain of small distances (Sz < d), the ng z dis- 
tributions are not Gaussian and their moments cannot 
be rescaled with a single exponent. In this range of Sz 
values, it would be possible to rescale all the moments 
{(h(z + Sz) — h(z)) p ) 1 l p by using several exponents £ p 
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FIG. 5: Probability density ng z (u) in a semi- logarithmic rep- 



resentation for 75 /im < Sz < 1.3 mm. Here u = with 
I = 35 /im and £ = 0.43 as obtained from the analysis of the 
second moment of the distribution (See Sec lIII A)l . For clar- 
ity only positive u are displayed. Continuous line : Gaussian 

1 -=° 2 
distribution —i=e 2 . 
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(not shown here): The relevance of such a multiscaling 
description at length scales for which the geometry of a 
single grain has a predominant effect is however question- 
able. 

These results bring new insight to a recently debated 
question concerning crack lines resulting from the rupture 
of p aper sheets that were found to exhibit multiscaling 
[23| . at least at small scales, i.e. at length scales compa- 
rable with the length of the fibers of the paper. Other 
observations, performed on the same material, but over 
a wider range of length scales, reported that, at larger 
scales, there was a crossover from this multiscaling be- 
havior towards a self-affine one 19]. These results may 
have similarities with our experiments, even though they 
correspond to the 3D rupture of a very different material: 
In our case, the crossover length correspond to the mean 
grain diameter d of the sandstone specimen. 



C. Dependence of the statistical properties on the 
sample width 

In addition to the lower limit of the order of the grain 
diameter d discussed above, the self-affine domain has 
also an upper limit £ (see Figure which generally de- 
pends on the sample size. For granite, the experimental 
upper limit £ is of the order of one fourth of the specimen 
size [H], HE)]. For mortar, the geometry can be consid- 
ered as self-affine only up to length scales of the order 
of 15 % of the specimen width W 18]. A systematic 
study of wood samples with various values of W leads to 
i = 0.1 W 

Figure [5] displays the variations of the correlation func- 
tion Ah(5z) as a function of 5z for samples of different 
widths W. All data sets display a self-affine domain fol- 
lowed by a plateau region corresponding to a value of 
Ah(Sz > £) increasing with W; the upper limit £ of the 
sclf-affinc domain also clearly increases with the system 
size. The variation of £ with W is plotted in Figure [7| 
and is well fitted by a straight line of slope 0.15 going 
through the origin. 

A value of the upper limit of the self-affine domain for 
profiles parallel to the direction x of crack propagation is 
also plotted in Figure [7] (0). This length is, in this case, 
the average of values of £ determined for four samples of 
different widths and the error bar reflects the variability 
of the values. In order to plot this result in Figure[7| these 
profiles have been assumed to correspond to a value of W 
equal to 35 mm {i.e. the length of the samples along the 
x direction). The corresponding data point is consistent 
with the linear fit corresponding in Figure [7] to profiles 
parallel to the z-axis. The scaling relation between £ and 
the system width W established above seems apparently 
also applicable in the x direction. 

Let us now compare the surface roughness in the self- 
affine domain at length scales between d and £ with the 
anomalous scaling characteristics recently observed 0, 
[I?], Oil]. Anomalous growth implies that, after a transient 
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FIG. 6: Log-log representation of the height— height correla- 
tion function Ah(Sz) computed along the z-axis for 4 spec- 
imens of different widths W. The straight line is a linear 
regression of slope £ = 0.47. 
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FIG. 7: Variations of the correlation length £ z (A) (resp. £ x 
(<>)) as a function of the specimen width (resp. length) W. 
Slope of the dashed line: 0.15. 



regime, the roughness amplitude Ah(Sz = est < £) scales 
as W^iobai-C w ith the system width W (see Eq.©). All 
the profiles analyzed correspond to the stationary prop- 
agation domain in order to make the comparisons mean- 
ingful. The extent of the transient domain near the crack 
initiation has been found in various materials to be of the 
order of few millimeters [H, [Hj]. Here, a value close to 
~ 2 mm was measured and all profiles considered in the 
present work are outside this region. 
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FIG. 8: Height— height correlation function Ah(8z) for fixed 
values of Sz as a function of the specimen width W. (o): 8z = 
1 mm; (o): 8z = 100 /im. Error bars: Standard deviation of 
Ah(8z). 



The height— height correlation function is plotted in 
Fig. [5] as a function of the specimen width W for two dif- 
ferent window sizes 5z. In contrast with the predictions of 
Eq. ([3]) for anomalous scaling, the roughness amplitude 
Ah(Sz < £) in the self-affine domain is independent of 
W over a range of variations of almost 1 to 10 (the same 
conclusions may already be inferred from Fig. [5] but the 
sensitivity is better here due to the linear scale). These 
characteristics differ therefore completely from those ex- 
pected for anomalous scaling. 



IV. DISCUSSION AND CONCLUSIONS 

Let us now discuss the origin of the statistical prop- 
erties observed for the present Fontainebleau sandstone 
surfaces. In a recent study, Bonamy et al. [l3j analyzed 
the fracture surfaces morphology in glassy materials, i.e. 
homogeneous silica glass and sintered glass. They suggest 
that, on the one hand, higher exponents £ ~ 0.8 reflect 
damage failure processes occuring in a zone surrounding 
the crack tips and are observable at length scales smaller 
than the size of the process zone; on the other hand, low 
roughness exponents £ ~ 0.4 would result from brittle 
fracture and be observable either in experiments where 
no process zone develops or at length scales larger than 
the size of the process zone. In this perspective, and as 
is discussed in more detail below, the low value Q ~ 0.45 
reported here for Fontainebleau sandstone samples is the 
signature of brittle fracture; moreover, since it is observed 
at length scales down to the grain size, no process zone 
is present in these experiments. 

In spite of these important differences, surfaces ob- 
tained from brittle and quasi-brittle failures share com- 



mon features. In both cases, a single exponent is enough 
to describe the scaling invariance of the roughness [19( : 
this confirms the relevance of the self-affine description 
of fracture surfaces for both types of materials. 

The lack of anomalous scaling in the present fracture 
surfaces compared to materials such as concrete, wood or 
granite displaying damage failure is an other indication 
that no damage process occurs in our experiments. In 
these latter studies, confinment effects limit the devel- 
opment of the process zone: Its extension and internal 
structure are then set by the system size W when a sta- 
tionary propagation regime is reached. As a result, from 
Eq. ([3]) [16J, the roughness amplitude Ah(Sz = est) varies 
with W as Ah{8z = est) oc H/ C gi°t>ai-C a s is indeed ob- 
served experimentally [2(|. In contrast, in the absence 
of a process zone, Ah(Sz = est) is only set by the mi- 
crostructural and mechanical properties of the material 
and does not depend on the system size (see for instance 
the model of Refs. [HI, [26]) as observed experimentally 
here. 

The same type of approach accounts for the domains 
of observation of the different self-affine geometries. For 
damage fractures, the upper bound £ of the self-affine 
domain where £ ~ 0.8 is set by the size of the process 
zone [l3, EH or by the sample size W if it is smaller 
than the process zone. In the present case of a brittle 
fracture with no process zone, the lower boundary of the 
self-affine domain is set by the grain size. The upper cut- 
off length £ increases linearly with the sample size W: in 
this case, no other characteristic length than W seems 
to determine £ so that it might increase indefinitely with 
W. These results suggest that the two different self affine 
geometries might coexist on a same fracture surface: the 
characteristic exponent would have a value C — 0-8 for 
lengths scales up to the size of the process zone and a 
value £ ~ 0.4 at larger length scales up to the plateau 
domain where Ah(8z) becomes constant. 

A key issue is whether the development of a pro- 
cess zone is purely material dependent. Recent exper- 
iments [13, suggest that fractured samples made of 
cemented grains, and therefore similar to Fontainebleau 
sandstone, do not always display brittle failure. The 
emission of acoustic waves was measured during the 
growth of a crack initiated from a chevron shaped notch 
under fatigue in a Permian sandstone sample. The lo- 
calization of the acoustic events during failure allowed 
to demonstrate that a process zone was present and to 
determine its size £ pz ~ 25 mm and its spatio-temporal 
evolution [27j • Moreover, the morphology of the result- 
ing fracture surfaces (28| was found to be self-affine for 
length scales less than 10 mm(< £ pz ) and characterized 
by an exponent ( — 0.75 ±0.03, similar to the usual value 
for materials displaying damage fractures. This observa- 
tion supports the assertion of Ref. [HI] and shows that 
the roughness exponent depends not only on the mate- 
rial, but also on the experimental conditions (geometry 
of initiation notch, loading mode...) 

To conclude, the present experimental results and their 
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comparison to previous works indicate that the value of 
the self-afhne exponent £ is a clear signature of the fail- 
ure mode : damage fracture involving a process zone if 
£ ~ 0.8 or brittle fracture if C, — 0.4. Another character- 
istic signature of damage fracture is anomalous scaling 
marked by a variation of the roughness amplitude with 
the sample size (if the latter is smaller than the process 
zone). The range of length scales over which the fracture 
geometry is self-afhne with a given exponent ( may rep- 
resent an important additional information. For £ ~ 0.8, 
this range should extend from the size of the grains to 
that of the process zone (which may be limited by the 
sample size); for £ ~ 0.4, it may extend from the size 
of the grains (if there is no process zone) or that of the 
process zone up to a length scaling as the size of the 
specimen. 

An important test of these hypothesis will then be to 
force the development of a process zone in materials of 
the type used here by modifying the loading or the ge- 
ometry of the experiments. For instance, one may use an 
initial chevron notch geometry as in refs. [U [Hj , where 
damage processes occur (so that C = 0.8) although the 



material is similar to the Fontainbleau sandstone used 
here. If a transition from brittle to damage fracture is 
induced in this way, this would result in a variation of 
the self affine exponent C or m the coexistence of two 
sclf-affinc domains with ( ~ 0.8 below the size of the 
process zone and £ ~ 0.4 above it. 
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